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Abstract. The formation of lithium lines in the atmo- 
sphere of C-rich giants is discussed. LTE and NLTE ap- 
proximations are used to model lithium lines in the spectra 
of super Li-rich AGB stars. The system of equations of the 
statistical balance of neutral Li in plane-parallel model at- 
mospheres is solved for a 20-level atom model. JOLA and 
line-by-line models of molecular absorption are used to 
compute synthetic spectra as well as the opacity in the 
frequencies of bound-bound and bound-free transitions of 
the lithium lines. Curves of growth and synthetic spectra 
are computed in LTE and NLTE for several model atmo- 
spheres of different Toff and C/0 ratios for four lithium 
absorptions, namely: the A4603, A6104, A6708 and A8126 
A Li I lines. The sensitivity of NLTE effects on Teg and 
the C/0 ratio is discussed. We found that NLTE mainly 
affects the resonance line doublet (A6708) while the im- 
pact of NLTE effects on the lithium subordinate lines, 
formed in the inner regions of C-rich giant atmospheres, 
is rather weak. Therefore the use of these lines is recom- 
mended for Li determination in AGB stars. However, in no 
case can we properly account for the formation of Li lines 
in AGB stars until sphericity, velocity stratifications, dust, 
chromospheres and other related phenomena, which are in 
fact present in AGB star atmospheres, are considered. Our 
results are used to derive Li abundances in three super- 
Li rich C-stars taking into account NLTE effects. Finally, 
the net Li yield from this class of stars into the interstellar 
medium is reconsidered. 
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1. Introduction 

The element lithium has provided many clues concern- 
ing stellar evolution and nucleosynthesis, as well as cos- 



mology. However, its origin is still far from clear. Nowa- 
days there is an increasing belief that this element has a 
multi-source nature: primordial nucleosynthesis (see the 
review by Wallerstein et al. 1997), galactic cosmic ray 
spallation in the interstellar medium and/or in the neigh- 
bourhood of supernova remnants (Meneguzzi et al. 1971; 
Feltzing & Gustafsson 1994), late type stars and novae 
(Abia et al. 1993a; Hernanz et al. 1996), supernovae ex- 
plosions (Woosley & Weaver 1995) and even spallation 
reactions around compact objects (Guessoum & Kazanas 
1998). However, except perhaps for the Li production dur- 
ing the Big Bang there is no consensus about the contri- 
bution of these sources to the present cosmic abundance, 
log e(Li)« 3.30 

The discovery of an unusually strong Li line at A6708 
A by McKellar (1940) in the carbon star (C-star) WZ Gas 
strongly suggested that AGB stars might be an impor- 
tant source of Li in the galaxy. The work of McKellar was 
followed by similar discoveries in other AGB stars: the 
C-stars WX Cyg and lY Hya (Sanford 1950; Abia et al. 
1991), the S-stars T Ara and T Sgr (Feast 1974; Boesgaard 
1970) and the SC-stars Henize 166 and VX Aql (Catch- 
pole & Feast 1971; Warner & Dean 1970). The same fig- 
ure has also been reported in AGB stars of the Magellanic 
Clouds (Plez et al. 1993; Smith et al. 1995). The pioneer- 
ing analysis by Cohen (1974) showed that the measured 
equivalent width of the A6708 A Li I line (Wa - 1 - 10 
A) in these stars is mainly the consequence of a strong en- 
hancement of the lithium abundance in the atmosphere. 
Quantitative determinations by spectral synthesis (Denn 
et al. 1991; Abia et al. 1993b; Plez et al. 1993; Smith et al. 
1995) have shown that the abundance of Li in these stars 
is 1-2 orders of magnitude higher than the present cosmic 
Li abundance. Therefore, these stars have been named su- 
per lithium-rich (SLiR) stars and might well constitute 
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^ The abundance of a given element X is noted as log e(X) = 
12-|-log(N(X)/N(H)) where N(X)/N(H) is the abundance by 
number of the element X. 
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the main source of Li in the galaxy. The production of 
Li in AGB stars has received special attention in recent 
years (e.g. Sackmann & Boothroyd 1992). Basically, Li 
is produced by the reaction '^Hc('*He, 7)^Be, followed by 
in a hot convective region that brings the 
^Be or ^Li to cooler regions before the ^Li is destroyed by 
'^Li(p, Q!)'*Hc. This is the so-called ^Be-transport mecha- 
nism (Cameron & Fowler 1971). Results of these theoret- 
ical studies quantitatively agree with the Li abundances 
derived in AGB stars of the galaxy and the Magellanic 
Clouds, although there are still many open questions re- 
lated to the Li production in these stars. Among others 
we find questions about the minimum initial stellar mass 
which can eventually become a SLiR star during the AGB 
phase, the duration of the SLiR phase or the actual Li 
yield into the interstellar medium. 

Our purpose in this work is related to the reality of the 
Li abundances in SLiR stars and its consequences on the 
Li yield. At present, the uncertainty in the derivation of 
the Li abundances in AGB stars is not lower than 0.4-0.5 
dex (see references above), mainly due to uncertainties in 
the stellar parameters and the fit to the observed spectra. 
Note, that in many situations no single set of stellar pa- 
rameters is found to fit the observations. However, there 
are also systematic errors that arc not visually taken into 
account that might dramatically change the Li abundance 
derived: uncertainties in the atmosphere models, the ex- 
istence of velocity stratifications (most SLiR stars arc ac- 
tually variable), sphericity and NLTE effects are not cur- 
rently considered as possible systematic sources of error. 
Of course, the consequences of each of these phenomena on 
Li abundance merit individual study and are beyond the 
scope of this work (see, however, Scholtz 1992; J0rgensen 
et al. 1992). Here, we will focus our attention on the ef- 
fects of departures from LTE in the formation of the Li 
lines in C-stars. 

The study of NLTE effects in the formation of lithium 
lines in stellar atmospheres was begun by the work of 
Miiller et al. (1975). They performed NLTE analysis of the 
very weak lithium resonance line A6708 A in the spectrum 
of the Sun. Later, Luck (1977) investigated the statistical 
balance of lithium in the atmospheres of G-K giants for 
a 4-lcvel atom model of Li. A similar atom model was 
used by de la Reza & Querci (1978) and de la Reza et al. 
(1981). Historically, these papers considered for the first 
time the impact of a stellar chromosphere on the lithium 
lines. A new stage of research was begun with the work 
of Steenbock & Holweger (1984), who used the technique 
of complete linearization for an 8-lcvcl atom model. They 
studied NLTE effects in lithium lines in atmospheres of 
dwarfs and giants. Pavlenko (1991) considered in detail 
the effects of deviation from LTE in the atmosphere of red 
giants. Later, Magazzii et al. (1992), Martin et al. (1994) 
and Pavlenko (1994), continued the NLTE studies in T- 
Tau stars, G-K giants, subgiants and dwarfs. The main 
results of these studies were confirmed by the indepen- 



dent work of Carlsson et al. (1994). who made a similar 
investigation using atmosphere models with various effec- 
tive temperatures, luminosities and metallicities. Finally, 
Houdebine et al. (1995), Pavlenko et al. (1995), Pavlenko 
& Magazzii (1996) and Martin et al. (1997) have contin- 
ued the studies of different aspects of NLTE formation of 
lithium lines in stellar atmospheres. 

As far as we know, the sole study on this subject for 
AGB stars is that by de la Reza & Querci (1978) who per- 
formed kinetic equilibrium calculations of neutral lithium 
lines in C-stars and determined the influence of the possi- 
ble chromospheric radiation into the photosphere. In the 
present work we revise this study using up to date atomic 
data, coUisional and radiative rates, an extended Li atom 
model and more reliable model atmospheres for C-stars in 
a wider range of effective temperatures (Toff = 2500—3100 
K) and C/0 ratios (1.0-1.35) (see below). We explicitly 
apply our results deriving Li abundances from synthetic 
spectra, both in LTE and NLTE, in three well known SLiR 
stars (WX Cyg, WZ Gas and lY Hya) from four accessible 
Li I lines: the resonance line at A6708 A and the subor- 
dinate transitions at A4603, A6104 and A8126 A, respec- 
tively, benefiting from the high signal-to-noise ratio and 
high resolution spectra of these stars. The consequences 
on the real Li abundances in AGB stars and on their net 
Li yield into the interstellar medium is then reexamined. 

2. Observations 

The observations were made during 1997 and 1998 in 

two different observatories. We used the 4.2 m WHT at 
the Observatory of El Roque de los Muchachos with the 
Utrecht Echelle Spectrograph as the main instrument and 
a 2048x2048 CCD with 24 fiia pixel size. We used the 
79.0 lines/mm grating which provides less wavelength cov- 
erage, but more space between orders (20-30 arcsec). The 
projected size of the slit on the chip was around two pix- 
els which gave a resolving power of 50000, the effective 
resolution ranging between 0.05-0.19 A from the blue or- 
ders to the red ones. The total number of orders on the 
chip were 30 covering the wavelength range 0.4-1.0 /xm 
with some gaps between orders. WZ Gas and WX Cyg 
were also observed by the 2.2 m telescope at the Calar 
Alto Observatory. For this observational run a fibre op- 
tics cassegrain echelle spectrograph (FOCES) (Pfeiffer et 
al. 1998) was used. This time the chip was a 1024x1024 
Tektronik CCD with 24 /xm pixel size. The FOCES image 
covers the visible spectral region from 0.38 to 0.96 /um in 
about 80 orders with full spectral coverage. Spectral or- 
ders are separated by 20 pixels in the blue and 10 in the 
red. The maximum resolving power is 40000 with a two 
pixel resolution element. 

The reduction of the spectra was made following the 
standard procedures using the ECHELLE task of the 
IRAF software package: bias subtraction, division by flat- 
field images, removal of the scattered-light, extraction of 



Abia et al.: Li formation in SLiR stars 
Table 1. Log of the observations and stellar parameters 



3 



St^T Spec, type Var. type %s C/N/0 ^'C/"C 



WX Cyg C8,2eJ M 2990 8.93/7.99/8.92 4.8 

WZ Cas"^ C9,2J SRa 3140 8.932/7.99/8.92 5.0 

lYHya N ? 2500 8.93/7.99/8.92 26 

WZ Cas has also been classified as a SC-type star. 
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the orders and wavelength cahbration with Th-Ar lamps. 
Typically the rms in the residuals of the calibrations were 
better than 10 mA. Finally, the spectra were divided by 
the spectrum of a hot rapidly rotating star to remove tel- 
luric absorptions, although they are only important in the 
A8126 A range. The signal-to-noise ratio of the spectra 
vary along the wavelength, from S/N'^ 20 in the bluest 
orders to S/N~ 700 around A8000 A. Since our stars are 
relatively bright objects (V^ 7 — 11), this strong variation 
in the S/N ratio achieved is mainly due to the very low 
emissivity of the stars below A ~ 4500 A. This is particu- 
larly evident in lY Hya, which is an extremely red object. 
For this reason, its spectrum in the A4603 A range is of 
too low quality (S/N~ 30) and has not been considered. 

3. The abundance analysis 

3.1. Atomic and molecular line lists 

The line list described in Lcbrc ct al. (1999) completed 
with C2 lines has been used for the analysis of the A6708 A 
region. Lines from the ^^€2, ^^C^^C and ^^C2 Swan sys- 
tem (Kurucz 1998) were considered and from the Phillips 
red system were predicted as in de Laverny & Gustafsson 
(1998). The line lists for the A4603,A6104 and A8126 A 
spectral domains were built in a similar way for the C2 
and CN molecules, and atomic line data were found in 
VALD (Piskunov et al. 1995). The broadening by radi- 
ation and van der Waals damping were calculated as in 
de Laverny & Gustafsson (1998). Line data were adjusted 
as described in Lebre et al. (1998) by comparison with 
the solar spectrum using the Holweger & Miiller (1974) 
model atmosphere for the Sun with element abundances 
from Anders & Grevesse (1989). A good fit to the solar 
spectrum was obtained in the A6708 and A8126 A spectral 
domains in ~ 50 A around the corresponding Li line. The 
fit, however, was not as good in the other two spectral 
ranges. Indeed, in the A4603 and A6104 A regions there 
are a number of unidentified features in the solar spec- 
trum with a modest absorption (~ 5 — 25 mA). Most of 
these features could be atomic in nature because no such 
intense molecular absorptions are expected in the Sun in 
these spectral ranges. However, we did not find any atomic 
lines in the data base of VALD or of Kurucz (1998) at the 
wavelengths of the missing features. Nevertheless, we were 
able to obtain a good fit to the solar spectrum in at least 
10 A around the A4603 A and A6104 A Li I features. Thus, 
we do not believe this problem will introduce an additional 
source of error into the abundance of Li derived from these 
two spectral domains. 

3.2. Atmospheric parameters 

Stellar parameters for the stars studied were taken from 
the literature when available. For WX Cyg the effective 

temperature derived by Ohnaka & Tsuji (1996) from the 
infrared flux method (IRFM) was adopted. For WZ Gas, 



effective temperature derived from the IRFM and angu- 
lar radii measurements (Dyck et al. 1996) agree quite well 
(±30 K). However, the value derived from these two meth- 
ods (3150 K) contrasts with that derived from infrared 
photometry (~ 2800 K; see Noguchi et al. 1981; Frogel et 
al. 1972). We adopted here the value obtained by Dyck et 
al. (1996). For lY Hya there is no estimate of its effective 
temperature nor its photometry is available in the litera- 
ture. We estimated the effective temperature by compar- 
ing its spectrum with the temperature sequence spectral 
atlas for C-stars created by Barnbaum et al. (1996). From 
that comparison we believe that lY Hya is a cool N-star 
and consequently, we adopted the lower effective temper- 
ature in our grid of atmosphere models i.e.: Tgff = 2500 
K. The uncertainty in Tefj is, however, no less than ±200 
K; it might well be larger for lY Hya. Note also that these 
stars are variable and we suspect Tgff variations as large 
as ~ 300 K during their cycle (see Richichi et al. 1995). 
The consideration of an unique effective temperature for 
these stars is therefore a strong assumption. 

A solar metallicity ([Fe/H]=0.0) was considered for the 
three stars. Most galactic C-stars are of Population I with 
near solar metallicity although typically they show strong 
enhancements of heavy elements (Zr,Ba,La etc...)(Utsumi 
1985; Dominy 1985). We checked that there is no intense 
heavy element line close enough to the four Li absorptions 
in our line lists that may affect the Li abundance deriva- 
tion. A study of the heavy element enhancements in the 
SLiR stars and other normal G-stars will be presented in 
a separate work. 

A gravity of log g= 0.0 and a microturbulcncc param- 
eter of ^ = 2.5 kms~^ were adopted. These are typical 
values for C-stars (see Lambert et al. 1986). Concerning 
gravity, note that studies by Pavlenko (1990) and Magazzii 
et al. (1992) on 0-rich dwarfs, subgiants and giant stars 
showed no important sensitivity in the ratio of LTE/NLTE 
Li abundances caused by changes in the gravity of Alog 
g= ±1, especially in the case of strong lines. Thus, uncer- 
tainties in gravity play a minor role in the formation of 
the lithium lines in G-rich atmospheres. Radiative damp- 
ing (which does not depend on log g) plays the main role in 
the formation of the wings of the saturated lithium lines. 
More serious seem to be the consequences of changes in the 
temperature structure of the atmosphere with decreasing 
gravity due to the increase in the effectiveness of spheric- 
ity effects. However, analysis of these effects lies beyond 
the scope of this paper. 

Finally, a model atmosphere was interpolated in Teg 
and the C/0 ratio for each star from an impublishcd grid 
of models for C-stars (Eriksson et al., private communica- 
tion). These models are constructed under the basic hy- 
potheses of hydrostatic equilibrium and plane-parallel ap- 
proximation. The models include more complete opacity 
data of HGN and C2H2 (in addition to diatomic molecules) 
with the sampling treated by an opacity distribution func- 
tion (see Eriksson et al. 1984, for details). The final C/0 
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Fig. 1. Theoretical LTE (solid) and NLTE (dotted) curves of growth for the four Li I absorptions for different Teff. 
The pairs of curves (solid-dotted) from top to bottom in each plot are for Teff=2500, 2800 and 3000 K, respectively. 



ratio of the model atmosphere used in the spectral synthe- 
sis was obtained by an iterative procedure comparing syn- 
thetic and observed spectra until a good fit was obtained. 
Li lines were not included in this procedure. However, the 
same C/0 ratio was not found to give the best fit in the 
four spectral domains for a given star, although differ- 
ences in the C/0 ratio never exceeded a few hundredths. 
For instance, in the A4603 and A6104 A regions we system- 
atically derived a lower C/0 ratio by ~ 0.03 hundredths. 



Table 1 shows the log of the observations and the final 
stellar parameters used in the analysis. The quoted CNO 
abundances are the mean values of those obtained from 
the best fit for each spectral domain (N abundance plays 
a minor role). 
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Fig. 2. Theoretical LTE (solid) and NLTE (dotted) curves of growth for the four Li I absorptions for different C/0 
ratios. The pairs of curves (solid-dotted) from top to bottom in each plot are for C/O=1.007, 1.02, 1.1 and 1.35, 
respectively. 



3.3. LTE abundances 

Lithium abundances were derived by synthetic spectra 
generated by using a modified version of code 
written at the Uppsala Astronomical Observatory. 
The theoretical spectra were convolved with the cor- 
responding FWMH to match the instrumental profile in 
the four spectral domains following the procedure used in 



Abia et al. (1993b). The carbon isotopic ratios {^"^C/^^C) 
derived in Abia & Isern (1997) from spectral synthesis to 
the 12 CN and ^^CN hnes in the A7990 - 8030 A spectral 
region were adopted (see also Abia & Isern (1996) and de 
Laverny & Gustafsson (1998), for details on the analysis 
and the choice of the CN parameters, respectively). Ta- 
ble 1 shows the carbon isotopic ratio adopted in the stars. 
They are in fair agreement with those derived by Ohnaka 
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& Tsuji (1996) and Lambert ct al. (1986) in WX Cyg and 
WZ Cas. The total error due to stellar atmosphere param- 
eter uncertainties, continuum placement and the fit itself 
amount to ±0.4 — 0.5 dcx for Li (LTE) abimdanccs and 
±6 — 13 for the ^^C/^^C ratio (see references above). 

Figures 1 and 2 (solid lines) show the behaviour with 
Toff and the C/0 ratio of curves of growth calculated in 
LTE for the four Li absorptions. As can be clearly seen, 
subordinate lines are saturated for high Li abundances 
(log e(Li)> 3). For the resonance line, however, the Li ab- 
sorption increases more steeply with increasing Li abun- 
dance. For very high Li abundances (log e(Li)> 4), the 
subordinate lines arc not very sensitive to the effective 
temperature, the C/0 ratio or gravity (as mentioned be- 
fore). This is a consequence of their saturated behaviour. 
The opposite occurs with the resonance line. This might 
lead us to consider the resonance line the better tool for 
abundance determinations in C-stars. However, the reso- 
nance line forms in the outermost layers (sec below) where 
the atmosphere structure is quite uncertain in C-stars and 
the microturbulence may even become supersonic. We will 
also see that NLTE effects are much larger for this tran- 
sition. 

Figures 3 to 6 show synthetic fits in LTE to the differ- 
ent spectral domains in WZ Cas. Synthetic spectra agree 
quite well with the observations for the A6708 and A8126 
A regions but this is not the case for the A4603 and A6104 
A ones. Some of the discrepancy is certainly due to the 
missing features in our line lists as mentioned in 3.1. Mo- 
roever, m the A4603 and A6104 A regions we obtain sys- 
tematically larger residual fluxes than observed. This, by 
the way, means that the Li abundances derived from the 
A4603 and A6104 A Li lines are systematically lower (see 
Table 3) than those derived from the red Li lines. The 
same figure was found by de la Reza & da Silva (1995) 
studying the formation of Li lines in K-giant stars. The 
reason for this is unknown. Bad gf-values for C2 lines, 
which are numerous in these bluest spectral regions, might 
partially explain this problem. Note that in these two spec- 
tral domains most of the C2 lines were taken from the data 
compilation by Kurucz (1998), which is probably not very 
accurate. Note also that we were not able to check their 
gf values with the solar spectrum as mentioned above. 
The discrepancy between theoretical and observed fluxes 
is particularly evident in the A4603 A region. We con- 
firmed that most of the strong absorptions in this spec- 
tral range are due to atomic lines while molecular lines 
contribute most to the background absorption. However, 
no reasonable reduction of the metallic abundances in the 
star can solve this discrepancy. Therefore, we are led to be- 
lieve that an incorrect figure for continuous opacity might 
be the main cause of this. 

The problem of the missing opacity in the blue has 
been widely discussed in the literature. Recently there 
have been some new results in this sense. Yakovina & 
Pavlenko (1998) showed that to fit the head of the strong 



NH band ^'^11 — AT'^S^ at A336 nm in the spectrum of the 
Sun, one should increase the continuum absorption coef- 
ficient in the region by a factor Ki, 1.7 — 1.9. Recently, 
Bell & Balachandran (1998) increased ^ 1.6 to flt Be II 
lines at A3 13.0 nm in the solar spectrum. We tried to arti- 
ficially simulate this missing opacity by increasing the H~ 
opacity by a free factor (namely, increasing the electronic 
pressure by a factor ~ 3). Indeed, in this case theoretical 
residual fluxes are lower (except for the very strong (satu- 
rated) lines) and the agreement between observed and the- 
oretical spectra improve. We preferred, nevertheless, to be 
coherent in our analysis in the four spectral domains and 
did not artificially change our set of conventional opacity 
sources. Thus, we used the same model atmosphere for all 
the spectral ranges and decided, for that reason, not to 
consider Li abundances derived from the A4603 A range. 
We will have to wait for the next generation of atmosphere 
models for C-stars to study this problem better (Plez et 
al. 1999). Note that in principle sphericity effects should 
be more pronounced in the blue part of the spectrum be- 
cause here the sensitivity to the temperature structure is 
higher even for strong spectroscopic features formed in the 
outer atmosphere. 

On the other hand, this discrepancy may be caused (at 
least partially) by the impact of dust opacity on the star's 
spectrum. Indeed, the formation of dust in the atmosphere 
of AGB stars and their envelopes has been known for a 
long time (see references in Wallerstcin & Knapp 1998). 
Dust particles can contribute to total opacity via scatter- 
ing processes, i.e. with cross-sections ~ 1/A". In that case 
the contribution of the dust opacity would increase to the 
blue. Li-rich K-giants also show an red excess, probably 
caused by a dusty envelope (although of a different na- 
ture; do la Reza & da Silva 1995). One may propose other 
explanations to explain this effect. For instance, the pro- 
cedure used to place the continuum level in the observed 
spectrum still seems very subjective. To provide more bal- 
anced conclusions regarding this problem we suggest that 
a more detailed study of these regions should be carried 
out in several SLiR stars. 

3.4- NLTE procedure 

To carry out the NLTE analysis for a 20-level Li atom 
model, we followed the procedure described in Pavlenko 
(1994) and Pavlenko & Magazzii (1996). A few specific 
items were taken into account here: 

— ionization-dissociation equilibria were computed for 
the carbon-rich case, i.e. considering C-contained 
molecules formation: CO, CH, CN, C2, CS, HCN, H2C, 
H3C, H4C, HC2, H2C2 and HC3. 

— in the frequencies of bound-free transitions of Li I 
the opacity of diatomic molecules was computed in 
the framework of just-overlapping-line approximation 
(JOLA, see Table 2). 
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Table 2. Data for JOLA opacities used in this work 



Molecule 


System 


gf 


Al (nm) 


A2 (nm) 


System names 


CaO 


C^E-XiS 


0.080 


730.0 


923.0 




cs 




0.011 


240.0 


330.0 




C2 


e^ng-a^n^ 


0.007 


237.0 


328.5 


Fox-Herzberg 


C2 


Ain„-xis+ 


0.011 


672.0 


1549.0 


Phillips 


C2 


b^E- - a^n,, 


0.002 


1100.0 


2700.0 


Ballik-Ramsay 


C2 




0.033 


340.0 


785.0 


Swan 


CO 


BiE+-Ain 


0.021 


412.0 


662.0 


Angsrom 


CO 


c^s+-A^n 


0.022 


368.0 


571.0 


Herzberg 


CO 


Ain-xis+ 


0.12 


114.0 


280.0 




so 


A^n-x^E- 


0.21 


246.0 


380.0 




SiO 


E1S+-X1S+ 


0.330 


171.5 


200.0 




SiO 


Ain-xis+ 


0.16 


207.0 


330.0 




NO 


c^n^-x^n^ 


0.0150 


207.0 


275.0 


6 


NO 


B^n^-x^n^ 


0.0035 


200.0 


380.0 




NO 


A^s+.x^n, 


0.0020 


195.0 


340.0 


7 


MgO 


BIS+-XIS+ 


0.048 


454.0 


544.0 




CN 


A2n-x2s+ 


0.0011 


400.0 


5630.0 


Red 


CN 


B2S+-X2S+ 


0.036 


240.0 


600.0 


Blue 


AlO 


c^n^ -x2s+ 


0.040 


200.0 


400.0 




AlO 


B2S+-X2S+ 


0.039 


404.0 


580.0 





Al and A2 aro tho limits for tlio molocular band opacity computations. 
Frank-Condon factors are from Kuz'menko et al. (1984) and gf values 
from Kuznetzova et al. (1980). 
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Fig. 3. Synthetic spectral fit to tlie A4603 A spectral domain in WZ Cas for log e(Li)=3.0. Note the strong discrepancy 
between observed (dots) and synthetic spectra in LTE (solid) and NLTE (dashed) (LTE and NLTE fits to the Li line 
almost coincide). Most of this discrepancy is probably due to an incorrect continuous opacity in this spectral range. 
The Li abundances derived from this line were not considered. 



— in the frequencies of bound-bound linearized transi- 
tions of Li I the atomic absorption was included. Fur- 
thermore, in the frequencies of the four Li lines we 
considered the absorption of CN and C2 for a given 
i2C/i=^C ratio (see Table 1). 

The work most relevant to the present study was done 
by de la Reza & Querci (1978; thereafter RQ). Using a di- 



vergence flux method these authors carried out an exten- 
sive modelling of NLTE effects of lithium lines in C-stars. 
They used a 4 level lithium atom model, atmosphere mod- 
els by Querci et al. (1974) and Johnson (1974) with effec- 
tive temperatures of 3800, 3600 and 3500 K (warmer than 
the Teft values usually derived in C-N stars), gravity log 
g= 1.0, C/0= 1.3 and modelled three Li I lines, namely 




Fig. 4. Synthetic spectral fit to tiie A6104 A spectral domain of WZ Cas for log e(Li)=3.0 LTE (solid line), NLTE 
(dashed line), observed spectrum (dots). 



the A6104, A6708 and A8126 A lines. They also tried to 
model the impact of a chromosphere on the formation of 
lithium lines using the approach of radiative temperature 
(Trad) on the bound- free transitions of lithium. Several 
aspects in their NLTE calculation differ from the present 
study: 

1. The atom model. We used a 20- level atom model, which 
allowed us to consider the interlocking of lithium lines 



(transitions) more appropriately. The ionization equi- 
librium of lithium is formed by the whole system of the 
bound-free transitions. Hovewer, as noted by RQ (see 
also Pavlenko 1991), transitions from/on the second 
level play the main role. 
2. Bound-free transitions. The radiation field in bound- 
free transition computations is very important in the 
modelling of lithium lines without LTE. Indeed, the 
effectiveness of the overionization of lithium depends 




Fig. 5. As Figure 4 in the A6708 A spectral domain of WZ Cas for log e(Li)=5.0 



directly on the mean intensities of the radiation field 

in the blue part of the spectrum (Pavlcnko 1991). 
The approximation T^ad = Te ( electronic temperature) 
used by RQ does not produce any overionization of 
lithium. Only in the case of T^ad > Tc docs the ion- 
ization of lithium increase with respect to the LTE 
case. Note that T^ad > Te was considered by RQ only 
when a chromosphcric radiation field was included. As 
noted by these authors, this approximation seems to 
be rather crude for the radiation field. 



3. bound-bound transitions. Lithium lines were treated 
by RQ as singlets. Basically, radiation transfer in 
the frequencies of several multiplet lines should dif- 
fer from the case one single (strongest) line. Moreover, 
wc include molecular line absorption in the continuous 
(background) opacity. In principle, this should reduce 
the probability of photon losses from the atmosphere, 
i.e. it directly affects the processes of radiative trans- 
fer. On the other hand, the radiative transfer in the 
continuum differs substantially from the case of con- 
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Fig. 6. As Figure 4 in the A8126 A spectral domain of WZ Cas for log e(Li)=4.5 



tinuum+lines in an atmosphere with chromosphere, as al. 1995) described in Pavlenko (1999)^. The de- 



here the impact of the hot chromospheric layers on the 
photosphere might be weaker. 



3. 5. NLTE curves of growth, synthetic spectra and 
abundances 



pendence of the departure coefficients hi — Nji^Lrpg/N^rpg, 
where NJ^lte ^lte ^'"^ t'^^ NLTE and LTE popu- 
lations of the i level, upon the depth in a typical C-star 
atmosphere is shown in Fig. 7. The behaviour of the de- 
parture coefficients of lower and upper lithium levels differ 
substantially. The populations of the lower levels are re- 



Comparison between LTE abundances derived with 
NLTE computations were performed by using a this code and those obtained with the Uppsala's code 
modified version of the code WITA2 (Pavlenko et showed an excellent agreement (~ ±0.01 dex). 
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duced by the ovcrionization processes. On the contrary, 
upper levels are overpopulated relative to the LTE case 
because they are more closely linked with the continuum 
(see also Pavlcnko 1994, Carlsson ct al. 1994, Pavlcnko & 
Magazzu 1996). Indeed, for log e(Li)> 4.0 the departure 
coefficients of the first level become very low in the out- 
ermost layers (6i < 10""* at Pg < 0.01), where the NLTE 
core of the Li resonance line forms. 

On the other hand, the depth of formation of lithium 
resonance and subordinate lines is quite different. For 
large lithium abundances (log e(Li)> 4.0) the formation 
region for the resonance line shifts toward the outermost 
layers, where the temperature drops below 1000 K. In fact, 
we have to extrapolate the model atmosphere structure to 
log Pgas ~ —5 to take the line formation consistently into 
account. This produces additional problems in the com- 
putation of coUisional and radiative rates. Here we are 
approaching the interstellar medium regime, where our 
approximations for coUisional and radiative rates may not 
be valid. Moroever, in this regime the splitting of Li terms 
into sublevels may also become inportant although, since 
we are dealing with saturated lines, we believe that this 
effect is not critical. Lithium subordinate lines are formed 
in deeper layers. For instance, even for log e(Li)> 4.5 the 
core of the A6104 A doublet forms in a region with temper- 
ature ~ 2000 K. As a consequence, the NLTE formation 
regime of this Li line changes dramatically: the interlock- 
ing processes of radiative transitions (sec Magazzu ct al. 
1992; Carlsson et al. 1994) become more important than 
overionization. In this situation, line profiles depend on 
the behaviour of the source function which is affected by 
the whole set of radiative transitions. In brief, considering 
the formation of lines in the outermost layers we note: i) 
the opacity here drops sharply due to the low molecular 
and H~ density, ii) the electron density is low and its effect 
on collision rates is small and iii) the electron temperature 
is also low, so that the effectiveness of the overionization 
oc e(~('^=~'^'''"^'/'^'^) would be high. 

Figures 1 and 2 (dashed lines) also show theoretical 
NLTE ciirvcs of growth for the four Li lines. Departures 
from LTE are very important for the resonance A6708 A 
line. NLTE corrections for this line can amount to 0.6 dex! 
However, NLTE effects for the subordinate lines are weak 
(< 0.2 dex) and even decrease for strong Li absorptions 
(high Li abundances, log e(Li)> 4). We also see from these 
figures that NLTE effects decrease for decreasing Tcff and 
increasing C/0 ratio in the model atmosphere. When 
Teff drops the opacity of the carbon-contained molecules 
increases so that the difference between Tj-ad (of the ra- 
diation field in the bound-free frequencies of Li) and Tg 
decreases. The same holds when the C/0 ratio increases in 
the atmosphere. It is interesting to note that LTE/NLTE 
abundances converge for high Li abundances for the sub- 
ordinate lines whatever the T^g and/or the C/0 ratio in 
the atmosphere. This gives us a chance to use subordinate 
lines as a tool for deriving Li abundances in Li-rich AGB 



stars even in the LTE approach, despite their saturated 
nature. Figures 3 to 6 show synthetic fits in NLTE (only 
for Li lines; i.e. dashed lines) in our stars. The difference 
between the LTE/NLTE fits to the resonance line is re- 
markable. For the subordinate lines, LTE and NLTE fits 
differ only slightly. 

Table 3 shows the final NLTE Li abundances derived in 
our stars. NLTE corrections are always positive, ranging 
from 0.1 to 0.5 dex. As was shown by Magazzu et al. (1992) 
for the case of G-M stars with solar abimdances, the signif- 
icance of NLTE effects depends upon the line strength. Li 
lines of moderate intensity (Wa > 0.2 A) form in a region 
where S^(rNLTE « 1) < B,^(tlte ~ !)• The NLTE cores 
of the resonance doublet become stronger than the LTE 
ones and, as a result, NLTE abundance corrections Alog 
e(Li)=log e(Li)NLTE— log e(Li)LTE are negative. However, 
in the case of the strongest (saturated) Li resonance dou- 
blets the difference between LTE and NLTE cores can- 
not be large (because both approach to zero) then, NLTE 
corrections again become positive as in the case of weak 
lines. Note that there is a region in the atmosphere where 
Si, < holds for all except the A4603 A Hue (Fig. 8) 
in the case of intermediate strong lines (log e(Li)~ 3). 
In these regions the interlocking processes of Li radiative 
transitions play the main role. 

4. Discussion and conclusions 

We now compare our LTE/NLTE results with the liter- 
ature. Abia et al. (1991) derived Li abundances in the 
same stars using the A6708 Li line. The abundances ob- 
tained here are considerably lower. In WZ Cas and WX 
Cyg the differences can be ascribed to the different choice 
of Tofi- and the C/0 ratio in the model atmosphere but 
also to the more complete molecular line list used in this 
work, which makes the background absorption more in- 
tense in the A6708 A region. Note that Abia et al. (1991) 
did not include the C2 molecule (in any isotopic form) or 
the ^^CN molecule in their line list. In lY Hya, moroever, 
the spcctriim now analyzed in the A6708 A range shows 
a much less intense Li line than that from which the Li 
abundance was derived by these authors: Weyos ~ 2.5 A 
here against 7.5 A in Abia et al. (1991). A similar varia- 
tion in time might be present in the spectrum of WX Cyg 
analyzed here (4.5 vs. 6.5 A). Whether this variation is 
an effect of a varying Toff along the pulsational phase or 
just evidence of the time scale for Li formation/depletion 
in AGB stars is a question which merits further studies. 
Note that Asplund et al. (1999) found the same figure 
in a post-AGB star and they concluded that they were 
observing the characteristic time scale of the mechanism 
responsible for the lithium production. 

Qualitatively our NLTE results agree with RQ but 
quantitatively some discrepancies are found. For in- 
stance, inspection to their Fig. 7 shov^^s a NLTE 
correction by ~ 0.7 dex at log e(Li)= 4 for the A6708 
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A line. However, our computations give a correc- 
tion of ~ 0.4 dex (see Fig. 2). This can be easily ex- 
plained by differences in the NLTE procedure, model at- 
mospheres etc. Namely, in this paper wc fully computed 
the opacities in the bound-free and bound-bound lithium 
transitions. This allowed us to establish the dependence of 
NLTE in lithium lines on C/0 and the ^^C/i^C isotopic 
ratio. Moreover, the analysis of RQ was based on the com- 
parison between the LTE and NLTE curves of growths. 



Our LTE /NLTE spectral synthesis approach to determine 
Li abundances in the atmosphere of Li-rich AGB stars 
seems to be better. Only in this way is the lithium abun- 
dance determination possible, since the conventional defi- 
nition of the equivalent width becomes meaningless due to 
the severe blending of Li lines. In this sense our results for 
subordinate Li lines arc of great interest. Since wc midcr- 
estimate the opacity of the bound-free and bound-bound 
transitions of lithium we get an upper limit to the NLTE 
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Fig. 8. Si//Bj/ ratio for the four Li lines in tiie atmospiiere 3000/0.0/1.007 of a C-star for different lithium abundances. 
Curves correspond to the A4603 (solid), A6104 (dotted), A6708 (long dashed) and A8126 A (short dashed) Li lines, 
respectively. For each pair of curves the upper line is computed with log e(Li)=3.0 and the lower one with log e(Li)=4.5. 



effects in these lines, but comparatively, we obtain low 
NLTE corrections on Li abundances. Using a more com- 
plete opacity source list, NLTE effects should be even less 
pronounced. That would comprise a good chance to obtain 
accurate abundance determinations based on the subordi- 
nate Li lines in the frame of the LTE approach, which is 
simpler and therefore more useful. On the other hand, 
we Ccin only qualitatively compare our NLTE re- 



sults with those by Pavlenko & Magazzu (1996) 
in M (O-rich) stars due to the obvious differences 
in opacities, model atmospheres etc, although in 
both cases, strong Li resonance lines are formed 
in the outermost layers of the atmosphere. Their 
curves of growth (see their Fig. 5) show a similar 
behaviour than ours (Figs. 1-2): for the A6708 A 
line W^™ > W^fLTE for fog e(Li)> 2; i.e., in the case 
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of strong lines formed high enough in the atmo- 
sphere NLTE effects reduce the equivalent width 
due to the dominant role of the overionization. For 

the A6104 A Li line they found W^te <W^lte foj. 
log e(Li)= 4.0 (see their Fig. 6). However, we ob- 
tain the contrary. This can be explained since in 
our case (C-rich atmospheres) the overionization 
should increase in the line-forming regions due to 
the drop of the density and pressure. 

Abia ct al. (1993a) computed an empirical Li yield 
from a homogeneous sample of galactic C-stars. They de- 
rived Li abundances from the A6708 A line in ~ 220 C- 
stars and by considering the mass-loss rate estimates for 
the stars in their study, obtained a Li yield into the in- 
terstellar medium of ~ 2 x 10"^ M0pc~^Gyr~^. These 
authors showed that this yield is extremely dependent on 
the stellar mass-loss rate assumed and on the real Li abun- 
dance in the star. Since approximately 90% of this value 
is determined by the actual yield from the SLiR C-stars, 
it is straightforward to revise this figure on the basis of 
the more accurate Li abundances obtained in this work. 
Now, considering the NLTE abundances obtained from 
the A8126 A Li line as the best estimate of the lithium 
abundance, we obtain the revised yield of ~ 8 x 10~^° 
Mopc~^Gyr^^, a factor ^ 3 lower than the previous value. 
This is mainly because the revised Li abundances in WX 
Cyg and lY Hya are considerably lower. As in Abia et al. 
(1993a), we can roughly estimate the contribution of C- 
stars to the galactic Li. Assume that the above production 
rate has been constant during the lifetime of the galaxy 
(~ 13 Gyr) and that the surface density of C-stars has 
been uniform and also constant during this time within 
a galactocentric radius of ^ 15 kpc (~ 50 C-stars kpc~^, 
see Claussen et al. 1987). In this case, the total contribu- 
tion to galactic Li abundance by C-stars is Mli ~ 3 Mq, 
i.e. ~ 10% of the total Li in the galaxy. This contribution 
might be higher if evolutionary effects are taken into ac- 
count, since the star formation rate was certainly higher 
in the past. Correspondingly, the number of C-stars that 
eventually became a SLiR C-star at a given time was also 
higher (see Abia et al 1993a for details). 

In summary: 

1. The formation of the lithium resonance line in C- 
stars is severely affected by NLTE effects. NLTE Li abun- 
dances are higher and can differ from the LTE ones up 
to 0.6 dex. The core of this strong Li line is formed in 
the outermost layers where we approach the interstellar 
medium regime. Furthermore, these layers are affected by 
different phenomena such as the stellar chromosphere, in- 
homogeneities, dusty shells, outflow/infall velocities etc. 
Our computations also show that lithium is severely ove- 
rionized. Therefore, the use of this saturated lithium line 
for abundance determinations in AGB stars seems quite 
impossible. In that sense, we confirm the results of RQ. 

2. Subordinate lithium lines are formed in the inner 
parts of the atmosphere. We have shown that NLTE ef- 



fects are rather weak even for strong lines (0.1-0.3 dex). 
Thus, they may be used for lithium abundance determi- 
nations in AGB stars even in the framework of the LTE 
approach. Of the three subordinate lines, the A8126 A Li 
I line is probably the best one for Li abundance deter- 
minations because it is less blended, because NLTE ef- 
fects are weak, because the molecular and atomic line list 
appears to be complete and because continuous opacity 
seems rather well reproduced in this spectral domain. 

3. NLTE effects in the lithium lines show a complicated 
dependence on the input parameters Tgff and the C/0 ra- 
tio. We found that NLTE effects increase with increasing 
Toff and a decreasing C/0 ratio. For a given TcS, NLTE 
effects decrease with increasing Li abundance for the sub- 
ordinate lines but increase for the resonance line. Further- 
more, for high Li abundances (log e(Li)> 4) LTE/NLTE 
abundances from the subordinate lines are not very sen- 
sitive to Teff or the G/0 ratio. For the resonance line, 
however, the opposite happens. 

4. The possibility of using the lithium subordinate dou- 
blet at A4603 A was considered. Theoretical spectra, how- 
ever, show a bad agreement with the observations, proba- 
bly due to a wrong estimate of the continuous opacity in 
this spectral range and to an incomplete line list in this 
region. Furthermore, this lithium line forms a strong blend 
with iron lines. Due to the very weak NLTE corrections 
found for it, the use of this Li doublet would be promising 
in the case of metal deficient stars if using very high res- 
olution spectra. Since we are confident of the atomic and 
molecular line list used beyond ~ 6000 A, we suggest that 
the C2 lines of the Swan system observed in the blue part 
of the spectrum are the largest source of error. More ac- 
curate line lists for these transitions are therefore needed. 
Such an improvement, together with better atmosphere 
modelling, solving the problem of the missing continuous 
opacity in the blue part of the spectrum, should improve 
the agreement between observed and synthetic spectra in 
the future. 

5. We have used observational data of high quality to 
determine lithium abundances in three C-stars: WX Cyg, 
WZ Gas and lY Hya. We have shown, that they are indeed 
SLiR stars with Li abundances in the range log e(Li)w 3 
to 5. However, up to now we have still failed to accurately 
determine Li abundances in AGB stars, mainly due to 
uncertainties in modelling their atmosphere. This leads to 
an important uncertainty concerning the estimate of the 
Li yield by AGB stars. These stars might account for up to 
30% of the currently observed Li or merely be a secondary 
source of Li in the galaxy. 
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Star A6104 A A6708 A A8126 A 

LTE NLTE LTE NLTE LTE NLTE 

WX Cyg 3.2 3.5 3.7 4.0 4.2 4.4 

WZ Cas 3.1 3.3 4.5 5.0 4.6 4.8 

lY Hya 2.6 2.7 2.5 2.7 3.2 3.3 
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